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A Note  on  Convolving 
Refinable  Function  Vectors 

C.  Conti  and  K.  Jetter 


Abstract.  When  convolving  two  refinable  function  vectors  which  give 
rise  to  convergent  subdivision  schemes,  the  convolved  scheme  is  again  con- 
vergent. Moreover,  the  conditions  on  the  mask  symbols  which  characterize 
the  approximation  order  of  the  associated  shift  invariant  spaces  show  that 
the  order  of  the  convolved  space  is,  essentially,  the  sum  of  the  order  of  the 
two  spaces  originating  from  the  convolution  factors. 


§1.  Kronecker  Convolved  Function  Vectors 

Our  previous  paper  [3]  deals  with  special  subdivision  schemes  associated  with 
a shift  invariant  space  of  bivariate  spline  functions,  where  the  “generators”  of 
the  shift  invariant  space  are  produced  through  convolving  lower  order  splines 
of  small  support.  The  present  paper  gives  a more  detailed  and  more  systematic 
analysis  of  this  convolution  process.  In  this  way  it  is  possible  to  prove  that 
(i)  the  convolution  of  convergent  subdivision  schemes  yields  a scheme  which  is 
convergent  as  well,  and  (ii)  essentially,  the  approximation  power  of  a convolved 
shift  invariant  space  is  at  least  the  sum  of  the  approximation  powers  of  the 
two  convolution  factors. 

We  recall  that  a vector  $ = 4>2,  ■ ■ ■ , 4>n)T  of  (continuous,  compactly 

supported)  d- variate  functions  is  called  refinable,  if  it  satisfies  a refinement 
equation 

$ = Aa  *(2--a)  . (1.1) 

a£7Z,d 

Here,  the  refinement  matrix  mask  A = {Aa)  ae7Zd  is  a matrix  sequence  with 
each  ‘coefficient’  Aa  being  a real  (n  x n)-matrix.  We  allow  only  masks  of  finite 
support,  i.e.,  Aa  =0  except  for  finitely  many  a € 2Zd. 
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Given  another  refinable  function  vector  = (ipi,rp2,  ■ ■ ■ ,i’m)T  of  d- 
variate  functions,  satisfying  the  refinement  equation 

* = £ Bp  *(2  • -/?)  (1.2) 

/3e2Z‘! 


with  corresponding  matrix  mask  B = g)  p&JLd , which  is  a matrix  sequence  of 

(raxni)- matrices  Bp,  we  use  the  following  Kronecker  type  notion  of  convolving 
the  two  vectors: 


© _ <J>  * ijr  ~ 


/<h**\ 
02**  ) 


\0n**/ 


Here,  the  convolution  of  a scalar  function  fa  with  the  vector  function  SS1 
is  taken  componentwise.  This  operation  produces  a vector  function  © = 
(#1,  (?2,  • . . , 0mn)T  with  m X n components  of  type  )m+j  '•=  <t> i * i>j-  It  is 
not  too  hard  to  see  that  © is  refinable  again, 


0 = £ Cy  0(2  ■ -7)  , (1.3) 

■V£7Zd 


where  the  refinement  mask  C = (C7)  which  is  a matrix  sequence  of 
(nm  x nm)-matrices,  is  computed  as  follows: 

^ £ A,  ® B^a  , a G 77/  . (1.4) 

aZTL* 

Here,  the  symbol  ® denotes  the  Kronecker  product  of  matrices. 

For  the  definition  and  some  properties  of  this  Kronecker  product,  we 
refer  to  [4,  Section  4.2].  It  should  be  noted  at  least  that  neither  the  Kronecker 
product  of  matrices  nor  the  Kronecker  type  convolution  is  commutative. 


§2.  Convergence  of  the  Convolved  Subdivision  Scheme 

Let  us  also  recall  that  a refinable  function  vector  gives  rise  to  a vector-valued 
subdivision  scheme  as  follows:  In  the  situation  of  (1.1),  the  subdivision  operator 
associated  with  the  refinable  function  vector  = (01,02, . . . , 0n)T  is  defined  as 

Sa  : (£(Hd))n  ->  (t(7Ld))n  , 

(Sa  A)„  :=  £ Al_w  Ap,  a G TLd  , (2-1) 

p€7Zd 

where  t(7Ld)  denotes  the  linear  space  of  sequences  indexed  by  7Ld . The  com- 
plete (stationary)  subdivision  scheme  consists  in  the  iterates  of  Sa,  namely: 

For  a given  initial  vector  sequence  A G (£(7Zd))n 
Put  A(°)  :=  A and 

Compute  A(fc+1)  :=  Sa  A^k\  k = 0, 1, . . . 


(2.2) 
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Hence,  the  iterate  A ^ = (Ao)1  )ae2Z,d  has  components 

A£fc)  :=  £ (4-a^)T  a?\  “ G ^ , (2-3) 


f>£7Ld 

with  the  iterated  matrices  A (*)  = (A^)^ Zi  defined  by  A^)  :=  A and 

:=  22  A0~1]Aa- 20  , ae7Ld,  for  k > 1.  (2.4) 


f)£7Li 


Following  [2,  Section  2.4]  we  say  that  the  subdivision  scheme  converges 
for  A = (A1, . . . , An)T  € (i°° (7Ld))n  if  there  exists  a continuous  function  fA  : 
]Rd  — > IR  such  that 


k^o  f^¥)e~A(k] 


= 0 for  e = (1  1 • • • l)r. 


(2.5) 


Here,  ||  • ll^  denotes  the  sup-norm  of  the  vector  sequence  A = (A1, . . . , \n)T 
given  by 

IIaII  :=  max  HA^oo. 

II  1 1 oo  i=l,...,n 

The  symbol  is  short  for  the  scalar- valued  sequence  (/a  (2^))agz;<i-  ^ 

should  be  noted  that,  since  for  any  convergent  scheme  the  limit  function  fA 
is  given  by 

/a  = 22  A°  $(-  - «)» 

aez:2 

we  can  recover  the  components  </>;,  i = 1, . . . ,n,  as  follows:  choose  the  initial 
sequence  A = (A1, . . . , An)T  in  the  following  way  that  A’  is  a delta  sequence 
(i.e.,  X’a  = 6a  for  a € TLd)  while  all  other  sequences  AJ,  j ± i,  are  null 
sequences.  This  special  initial  vector  sequence  and  its  iterates  will  be  denoted 
by  Ej  and  E k > 0,  respectively.  It  then  turns  out  that 


^Jk%)e-E-1L=o  for  (2-e) 

Theorem  2.1.  Given  two  convergent  subdivision  schemes  associated  with  the 
refinable  function  vectors  $ and  4/,  the  (Kronecker  type)  convolved  scheme 
associated  with  © = * fir  is  again  convergent. 

Proof:  Let  f :=  (11  •••  1)T,  and  let  F j denote  the  vector  sequence  (com- 

771 

posed  by  m sequences)  with  the  delta  sequence  at  position  j and  the  null 
sequence  at  all  other  positions.  If  k > 0,  are  the  iterated  vectors  with 
respect  to  the  ^-subdivision,  we  have 


^(^)f-F f =0  for  j = 1 


. . ,m, 


oo 


(2.7) 
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in  addition  to  (2.6).  In  order  to  prove  the  theorem,  we  will  show  that 


lim 


(fell 


= 0 


(2.8) 


for  i = 1,  ...,ra  and  j = Here,  the  convolution  of  two  vector  se- 

quences is  defined  by 


(xl\ 

/71  \ 

/a1  *r\ 

A2 

* 

72 

a2*t 

Wn) 

\7m/ 

\ A"  * r / 

where  the  convolution  of  a scalar  sequence  A*  with  the  vector  sequence  T is 
taken  componentwise,  i.e.,  A*  * T = (A*  * 71,  A'  * y2, . . . , A1  * ym)T . We  will 
also  use  the  estimate 


l|A*r||00  < min  { ||A||i||r||oo  , ||A||  Jr^  }, 
where  ||A||  :=  max^^,..^  ||A'||i,  with  the  usual  1-norm  of  scalar  sequences. 

Now,  the  iterated  matrices  C ^ for  the  ©-subdivision  can  be  expressed 
by  the  iterated  matrices  A W and  for  the  $-  and  ^-schemes  as  follows: 

- 2^  E A*k)  ® B*-s  - “ e E"-  (2-9) 

Thus,  taking  E*  * F,-  as  a starting  vector  for  the  ©-subdivision,  the  iterated 
vectors  are  given  by 

(E»  * Fj)(fc)  = E<fc)  * Ff  , k > 0,  (2.10) 

whence 

((^*^K2fc))(e®f)-(Ei*P,-)(fc) 

= ((**^j)(2fc))(e®f)- ^(^(^e)  * 

+ 2SfEi‘)*(^(2F)f-F?)) 

We  estimate  the  three  vector  sequences  in  the  preceding  three  lines  as 
follows:  the  first  term  is  a vector  sequence  where  each  component  of  the 
vector  consists  of  the  sequence  (e^)agZd  with 

e“fc)  :=  L (^(^)) 

” SL/-r77d  ~ 
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This  is  the  error  of  a tensor  product  rectangular  rule  applied  to  the  convolution 
integral,  and  due  to  the  continuity  and  the  compact  support  of  <pi  and  tpj  v/e 
get 

lim  sup  |e£^|  = 0. 
k—00aeHd 

For  the  third  term, 


1 

2 d'k 


and  the  bound  tends  to  zero  as  k — > oo,  since 

= a?*  z h(?)|  ^ ^ \Mx)\dx- 

Finally,  for  the  second  term 


ct£TLd 


rafE?’  * < « ||*<2.)f-E 


2d ■ 


as  k — > oo.  This  follows  from 


Ck  2dfc 


E 


(fc) 


< — r-r  0i(— r)e  H — — r)e  — 

l-  2d'k  y2k  ’ i 2d  k Viy2k>  1 l 


since  the  first  term  on  the  right-hand  side  converges  to  JRd  \<f>i{x)\dx,  and  the 
second  term  is  bounded  due  to  (2.6),  the  uniform  continuity  of  the  compactly 
supported  function  fa,  and  the  fact  that  for  some  compact  set  K c !Rd  (in 
fact  we  can  take  the  support  of 


(Ef  ))„  = 0 for  ~iK. 

This  latter  property  is  a consequence  of  the  compact  support  of  the  matrix 
mask  A,  since  we  start  the  iteration  with  a ‘delta’-sequence  E,.  In  conclusion, 
(2.8)  holds  true.  □ 


§3.  Approximation  Order 

Approximation  orders  of  shift  invariant  spaces  have  been  studied  quite  inten- 
sively in  the  past  few  years.  Concerning  definitions  and  notation  we  refer  to 
the  recent  survey  [5].  A characterization  of  approximation  power  in  terms  of 
the  mask  symbol  is  also  given  there.  In  case  of  the  refinement  equation  (1.1), 
this  mask  symbol  is  defined  by 

h*(0  :=  Y*  £ A° 

a£2Zd 


(3.1) 
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The  result  is  as  follows:  For  given  k 6 IN  we  say  that  H$  satisfies  condition 
( Zk ) if  there  exists  a row  vector  v = (tj  , . . . , rn)  of  trigonometric  polynomials 
such  that 

v(0)  4>(0)  ± 0,  v(0)H*(0)=v(0),  (3.2 .a) 

and 

=0  for  \p\  < k and  0 € E'0,  (3.2.6) 

where  E'0  is  short  for  the  set  of  corners  of  the  cube  [0,  l]d  with  the  origin 
removed.  [5,  Theorem  3.2.8.(i)]  then  asserts  that  condition  (Zk)  implies  that 
S$  has  I<2-  approximation  order  k,  for  f € W*  (Rd). 

Lemma  3.1.  If  the  mask  symbol  satisfies  condition  (Zk)  (with  the  row 
vector  v)  and  the  mask  symbol  H^,  satisfies  condition  (Z/)  (with  the  row  vec- 
tor w ) then  the  mask  symbol  H®  of  the  convolved  function  vector  © :=  <i>  * ’3> 
satisfies  condition  Zk+i  (with  the  row  vector  z :=  v ® w). 

Proof:  We  apply  the  identity 


H®  = H$  ® H*  (3.3) 

several  times.  Condition  (3. 2. a)  holds  for  the  convolved  function  vector,  since 
(v  ® w)(0)  (3>  * ^)(0)  = v(0)4>(0)w(0)4»(0)  ^ 0, 


and  by  (3.3), 

z(0)H©(0)  = (v(0)Hf(0))  ® (w(0)H*(0))  = v(0)  ® w(0)  = z(0). 

In  order  to  verify  condition  (3.2.6), 

jDm(z(-)H©(-))  |=0  for  \p\  <k  + i and  f3  € E'0, 

we  make  use  of  the  Leibniz-type  formula 

^((v®  w)(-)(H*®H*)(-)) 

= E (!)  W (V()H* (-))  ® D (w(-)H,(l)) 
o<M<H  V7' 

at  every  point  (3  € E'0.  □ 

In  order  to  derive  an  approximation  order  result  from  this  lemma,  we  refer 
to  the  precise  statement  of  [5,  Theorem  3.2.8].  A sample  result  is  as  follows: 
If  the  shift  invariant  spaces  associated  with  $ and  \F  have  approximation 
orders  k and  l,  respectively,  and  if  the  Gramians  and  Gq,  satisfy  the 
regularity  condition  given  in  [5,  Theorem  3.2.8.(ii)],  then  the  ‘convolved’  shift 
invariant  space  has  (at  least)  approximation  order  k + l.  However,  in  general 
the  ‘convolved’  Gramian  G©  does  not  satisfy  this  regularity  condition. 
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§4.  An  Example:  Bivariate  C 1 Cubics  on  a 3-directional  Mesh 

In  [3],  we  have  given  an  example  of  piecewise  C^-cubics  on  the  four-directional 
mesh.  Following  [1],  we  consider  piecewise  C1-cubics  on  the  three-directional 
mesh  generated  by  the  lines  x = k,  y = l,  with  k,  l 6 2Z  when  adding  the 
diagonals  x—y  = m,  m £ TL.  A basis  of  this  space  is  given  by  the  two  functions 
0 1 = Bm *XTi  and  62  = Bm*XT2,  with  B m the  linear  three-directional  box- 
spline  (or  “Courant”  element),  and  xt2 , Xt2  the  characteristic  functions  of  the 
two  triangles  Ti,T2  obtained  by  cutting  the  unit  square  [0,  l]2  by  the  ‘north- 
east’ diagonal.  Thus,  © = $ * with  $ :=  (Bin)  and  :=  (xt1,Xt2)T  ■ 
Now,  $ = (Bm)  satisfies  a scalar  refinement  equation  (1.1)  with  refine- 
ment mask 


A — (Aa)ae7Z2  — 


0 0 0 0 0 


0 0 
0 1 


2 2 

1 i 


oil 

u 2 2 


0 

5 0 
0 0 


00000 


V 


J 


and  'V  = (xt,  , XTi  )7  satisfies  the  refinement  equation  (1.2)  with  matrix  mask 


/ 


B — (Ba)ag2 z2  — 


\ 


0 0 
0 0 

0 0 
0 0 


( 0 0 
VI  1 


1 0 
0 1 

0 


0 

1 0 
0 1 

1 1 
0 0 

0 


Here,  the  indexing  of  the  ‘coefficients’  is  such  that  the  boldface  entry  is  at 
position  a = (0, 0). 

It  follows  that  © = satisfies  the  refinement  equation  (1.3),  with  the 
matrix  mask  C as  displayed  on  the  following  page,  and  Theorem  2.1  provides 
the  convergence  property  for  the  associated  subdivision  scheme. 

Concerning  approximation  order,  Lemma  3.1  can  be  applied  by  putting 
k = 2,  i = 1,  and  z :=  (1)  <g>  (1  1).  As  a consequence,  the  approximation  order 
for  bivariate  C'1-cubics  on  the  the  three-directional  mesh  is  at  least  3,  and  this 
is  the  precise  approximation  order  (as  was  proved  in  [1]). 
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